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There  are  two  elements  in  this  project.  They  are, 

1.  Analysis  of  the  noise  of  high-performance  military  aircraft. 

2.  Investigation  of  the  importance  of  indirect  combustion  noise  as  a  dominant  component 

of  military  aircraft  noise  at  afterburner. 

This  quarterly  progress  report  concentrates  on  the  second  element  of  the  project. 

I.  Introduction 

When  extra  fuel  is  burnt  at  afterburner  operating  condition,  the  combustion 
process  is  highly  unsteady.  This  creates  large  temperature  fluctuations  resulting  in  the 
generation  of  numerous  high  and  low  temperature  blobs.  These  blobs  are  often  referred  to 
as  entropy  waves.  It  is  now  known  that  when  entropy  waves  or  entropy  blobs  are 
convected  through  a  mean  flow  with  significant  axial  velocity  gradients,  indirect 
combustion  noise  is  generated  (see  figure  1).  The  present  research  is  to  investigate  this 
phenomenon  computationally. 


Figure  1.  The  passage  of  hot  and  cold  spots  of  combustion  from  the  afterburner 
through  the  convergent-divergent  nozzle  is  the  mechanism  responsible  for  the 
generation  of  indirect  combustion  noise. 


Two  critical  elements  are  required  in  any  study  of  indirect  combustion  noise 
generation  by  numerical  simulation  in  a  military-styled  nozzle.  They  are, 

1 .  The  creation  of  a  random  broadband  entropy  wave  field  as  an  input  to  the  computation. 

2.  A  mean  flow  inside  the  nozzle. 
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Once  these  two  elements  are  established,  a  parametric  study  of  indirect  combustion  noise 
can  be  carried  out  by  passing  the  computed  broadband  entropy  wave  field  through  the 
mean  flow  inside  the  nozzle. 


In  this  progress  report,  we  will  report  the  following. 

A.  The  development  of  a  stochastic  entropy  wave  generation  model  to  be  used  as  input 
boundary  condition  in  a  numerical  simulation.  We  will  first  report  the  development  of 
a  one-dimensional  model  as  it  is  less  complicated  and  hence  easier  to  understand.  Our 
report  below  includes  the  analytical  model  and  validation  tests.  Then  we  report  the 
extension  to  three-dimensions  and  validation  test  results. 

In  developing  our  three  dimensional  stochastic  entropy  wave  generation  model,  we 
prescribe  the  requirements  that  a  user  must  be  able  to  specify, 

a.  The  intensity  of  temperature  fluctuations. 

b.  The  spectrum  of  fluctuations. 

c.  The  two-point  cross  correlation  functions  in  two  lateral  perpendicular  directions  in 

the  lateral  plane  perpendicular  to  the  direction  of  the  mean  flow.  These  two-point 
temperature  fluctuation  cross-correlation  functions  are  designed  to  allow  a  user  to 
specify  the  lateral  size  of  the  entropy  blobs. 

B.  The  development  of  a  high  quality  computational  algorithm  for  computing  the  mean 
flow  inside  a  military  styled  nozzle.  This  algorithm  permits  a  user  to  prescribe  all  the 
mean  flow  variable  at  the  inflow  boundary.  They  include  the  pressure,  temperature 
and  density  in  addition  to  the  inflow  velocity  and  Mach  number.  One  characteristic 
feature  of  a  military-styled  nozzle  is  that  it  has  a  fairly  sharp  throat.  Because  of  the 
sharp  throat  shocks  are,  invariably,  created  inside  the  nozzle.  Our  algorithm, 
therefore,  is  designed  with  shock  capturing  capability. 

II.  Stochastic  model  boundary  condition  simulating  incoming  broadband  entropy 
waves/blobs 

A.  One  Dimensional  Stochastic  Entropy  Wave  Field  Generation  Model 

We  will  use  dimensionless  variables  with  the  following  quantities  as  scales 
A  (the  size  of  the  dominant  entropy  blobs)  as  the  length  scale. 

aQ  (speed  of  sound)  as  the  velocity  scale 
A  /  aQ  as  the  time  scale 
po  as  the  density  scale 
p0a20  as  the  pressure  scale 

For  one  dimensional  compressible  flow,  the  governing  equations  for  entropy  and 
acoustic  waves  in  dimensionless  forms  are: 
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where  an  overbar  indicates  a  mean  flow  quantity. 

We  will  let  the  spectrum  of  the  prescribed  temperature  fluctuations,  5 (to) ,  be  as 
shown  in  figure  2. 


Figure  2.  Discretization  of  the  prescribed  spectrum  into  multiple  bands  with  width  A co. . 

We  propose  to  discritize  the  given  spectrum  S (to)  into  multiple  narrow  bands 
with  center  frequencies  co.  ( j=l  to  N)  and  bandwidth  A  co.  as  shown  in  figure  2.  The 
width  of  the  narrow  bands  Ato .  is  random. 


It  is  straightforward  to  show  that  the  following  stochastic  temperature  field 
satisfies  governing  equations  (1),  (2)  and  (3). 
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We  may  consider  equation  (4)  as  consists  of  multiple  simple  harmonic  oscillators  with 
frequencies  at  co.  .  To  conserve  energy,  the  amplitude  A  of  each  oscillator  with 

frequency  co.  is  given  by, 

±/r  =  2s(ffl)A®  (5) 


To  randomize  the  oscillators,  in  this  model,  the  phase  <j>  is  to  be  a  random  number  scaled 

to  the  range  of  [o  to  2/rJ .  So  the  random  wave  field  given  by  equation  (4)  is  created 
primarily  by  phase  randomization. 

To  show  that  the  model  given  by  equation  (4)  would  yield  a  stochastic  field  with  a 
spectrum  and  intensity  the  same  as  the  prescribed  input  i.e.  j  of  figure  2,  we  can  first 

compute  the  two-point  space-time  correlation  function  from  the  stochastic  field  generated 
by  equation  (4).  By  means  of  equation  (4),  after  some  lengthy  algebra  and  the  repeated 

use  of  the  formula  cosAcos5  =  ^cos(i4-5)+cos(A-i-5)J  to  reduce  the  product  terms 

of  cosines  to  a  sum  of  cosines,  it  is  straightforward  to  derive, 
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where  the  overbar  denotes  a  time  average.  Also, 
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where  8jk  is  the  Kronecker  delta.  On  substituting  (7)  into  (6),  the  space-time  cross¬ 
correlation  function  is  found  to  be, 


T(x',t')r(x",t'+T)=  £  s(tojA®fccos 
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The  autocorrelation  function  of  the  entropy  field  may  be  found  by  setting  x"  =  x' .  This 
yields, 


-  N 

r2(x',t')  =  ^  — >  J  S^co^dcQ-,  in  the  limit  Acok  — >  0  (9) 

k=-N  _o» 

Thus  the  spectrum  of  the  stochastic  entropy  wave  field  generated  by  equation  (4)  will 
yield  a  spectrum  equal  to  s(w  j ,  the  input  spectrum. 


Validation 

To  illustrate  the  effectiveness  of  our  model,  let’s  consider  the  following  test  case. 
Consider  a  one  dimensional  duct  carrying  a  uniform  mean  flow  of  Mach  number  M.  A 
stochastic  entropy  wave  field  enters  the  duct  on  the  left  end  and  exits  on  the  right  end.  Let 
the  size  of  the  dominant  hot  spot  be  1  in  dimensionless  unit.  We  would  like  to  be  able  to 
resolve  hot  spots  0.2  times  smaller  that  the  dominant  ones.  To  be  able  to  do  so,  we  will 
take, 

Ax  =  0.2/7 -0.02 

The  factor  1/7  above  comes  from  the  resolution  of  the  DRP  scheme  that  we  intend  to 
use  in  all  our  computations. 

The  entropy  blobs  are  convected  downstream  by  the  mean  flow.  So,  in  one  unit  of 
dimensionless  time,  the  flow  moves  a  distance  of  M  past  a  fixed  point.  The  number  of 
dominant  hot  spots  passing  this  point  per  unit  time  is  M  /1.0  =  M  .  Therefore,  the 
dominant  frequency  of  entropy  blobs  passing  a  point  in  the  computational  domain  is  M . 
Hence,  the  dominant  angular  frequency  is  coQ  =  2 nM  .  We  will  use  a  Gaussian  spectrum 
with  the  following  parameters. 
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Figure  3.  An  instantaneous  temperature  distribution  in  the  test  duct  taken  from  a 
numerical  simulation. 


To  demonstrate  the  capability  and  accuracy  of  the  model  a  computer  code  is 
developed  to  generate  the  stochastic  broadband  wave  field  according  to  equation  (4).  The 
computed  space-time  data  of  temperature  is  processed  for  comparison  with  model  input 
functions.  Figure  3  shows  a  computed  instantaneous  temperature  distribution  in  the 
computation  domain.  The  temperature  field  varies  in  space  and  time.  There  are 
approximately  10  blobs  in  the  computational  domain.  This  is  expected,  as  the  dominant 
entropy  blob  has  a  length  of  1.0  in  this  test  case.  Figure  4  is  a  plot  of  the  input  spectrum 
(full  line).  Also  shown  is  the  spectrum  computed  from  the  numerical  data  of  the 
simulation  (dashed  line)  measured  at  x  =  5.0.  Figure  5  shows  the  input  autocorrelation 
and  the  autocorrelation  obtained  by  processing  the  numerical  data.  It  is  self-evident  that 
there  are  good  agreements  between  the  input  and  numerical  results  in  both  figures  4  and 
5.  This  provides  support  to  the  belief  that  the  entropy  wave  generation  model  developed 
in  this  project  is  accurate  and  reliable. 
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Figure  4.  Comparison  between  input  and  computed  power  spectral  density  of  the 
temperature  fluctuations  of  the  test  example.  Input  (full  line),  processed  from  numerical 
data  (dashed  line) 


X 


Figure  5.  Comparison  between  input  and  computed  autocorrelation  function  of  the 
temperature  fluctuations  of  the  test  example.  Input  (full  line),  processed  from  numerical 
data  (dashed  line) 

Dr.  Schuster  of  Honeywell  provided  us  with  his  temperature  fluctuation  spectrum 
measured  near  the  exit  of  the  combustor  of  the  Honeywell  TECH977  turbofan  engine. 
The  spectrum  is  shown  as  a  dashed  line  (reddish  color)  in  figure  6.  The  smooth  curve  in 
this  figure  is  the  corresponding  spectrum  computed  from  data  generated  by  our  model 
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using  Dr.  Schuster’s  spectrum  as  input.  There  is  excellent  agreement  over  the  entire 
frequency  range  up  to  4000  Hz. 


Figure  6.  Temperature  fluctuation  spectrum  measured  by  Schuster  et  al  in  their  engine 
noise  study  using  TECH  977. 

B.  Three  Dimensional  Stochastic  Entropy  Wave  Field  Generation  Model 

We  have  successfully  generalized  the  one-dimensional  model  of  equation  (4)  to 
three  dimensions.  In  three  dimensions,  in  addition  to  prescribing  the  intensity  and 
spectrum  of  the  random  wave  field,  a  user  can  prescribe  the  two-point  cross  correlation 
functions  in  two  perpendicular  directions  in  the  lateral  plane.  The  half-widths  of  these 
two-point  cross  correlation  function  serve  as  a  measure  of  the  lateral  size  of  the  random 
entropy  blobs,  Thus,  if  the  x-axis  is  the  direction  of  the  mean  flow,  then  one  can  prescribe 
the  spectrum  and  the  two-point  cross  correlation  functions  in  the  y  and  z  directions. 


Figure  7.  Inclined  plane  at  30°  to  the  mean  flow  direction  for  observation  of  stochastic 
entropy  wave  field. 
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Figure  8.  Stochastic  entropy  wave  field  generated  by  3-D  model  with  spectrum  and 
lateral  cross  correlation  function  given  by  equations  (1 1)  to  (13). 

To  illustrate  the  accuracy  of  the  model,  we  consider  the  test  case  with  the 
following  spectrum  and  lateral  cross  correlation  functions. 


n  J 


2(7 


(11) 


T(x,y',z,t)T(x,y",z,t)  = 


y'-y" 

b 


(12) 


T(x,y ,z' ,t)T(x,y ,z" ,t)  =  T2e  '  b  ’ 
The  parameters  are  assigned  the  values  below. 

T2  =  1.0,  b  =  0.75,  M  -0.3,  coQ  -  2nM,  o  =  cqJ2 


(13) 


To  see  the  entropy  wave  field  in  3  dimensions,  we  consider  the  plane  OABC  in  figure  7. 
The  mean  flow  is  in  the  x-y  plane  at  45°  to  the  y-axis.  Figure  8  shows  the  entropy  wave 
field  generated  by  the  model  as  observed  on  the  inclined  plane  of  figure  7.  In  figure  8,  red 
indicates  high  temperature  and  blue  implies  low  temperature. 
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y 

Figure  9.  Comparison  between  input  two-point  cross  correlation  function  and  the 
computed  cross  correlation  function.  Here  y  -  [y'—y") . 


z 

Figure  10.  Comparison  between  input  two-point  cross  correlation  function  and  the 
computed  cross  correlation  function.  Here  z  =  [z'~  z") . 
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Figure  9  shows  a  comparison  of  the  two-point  cross  correlation  function  in  the  in¬ 
direction.  The  full  line  is  the  input  function  (equation  12).  The  black  dots  are  the  cross 
correlation  function  computed  from  the  model  generated  stochastic  temperature 
fluctuation  field.  There  is  excellent  agreement.  Figure  10  shows  similar  comparison  in  the 
z-direction.  Again  there  is  excellent  agreement.  Figure  1 1  shows  a  comparison  between 
the  input  spectrum  function  and  the  spectrum  function  computed  from  model  generated 
random  field  data.  The  agreement  is  again  excellent.  The  excellent  agreements  provide  a 
strong  validation  of  the  3  dimensional  model  boundary  condition  developed  in  the  present 
research  effort. 


Figure  11.  Comparison  between  input  temperature  fluctuation  spectrum  and  the 
spectrum  computed  from  stochastic  temperature  fluctuation  field  generated  by  the  model 
boundary  condition  developed  in  the  present  research  effort.  Red  line  =  input  spectrum. 
Black  dots  =  computed  spectrum. 


III.  Computation  of  mean  flow  inside  a  military  styled  nozzle 

Recently,  Tam  and  Parrish  ( Noise  of  High-Performance  Aircraft  at  Afterburner 
Journal  of  Sound  and  Vibration,  352,  103-128,  2015)  proposed  that  indirect  combustion 
noise  could  be  an  important  noise  component  of  the  F-22  aircraft  at  afterburner.  F-22 
aircraft  have  rectangular  nozzle.  To  ensure  a  high  chance  of  success,  we  will  first 
investigate  the  generation  of  indirect  combustion  noise  inside  a  rectangular  nozzle.  Figure 
12  shows  the  upper  half  of  the  nozzle  geometry  and  Mach  number  to  be  used  in  the 
present  work.  The  exit  Mach  number  1.4  is  based  on  a  one  dimensional  compressible 
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flow  calculation.  For  three  dimensional  flow,  shocks  will  appear  inside  the  nozzle  so  that 
the  flow  at  the  nozzle  exit  will  not  be  uniform. 


Figure  12.  Geometry  of  the  upper  half  of  the  rectangular  military  styled  nozzle  used  in 
the  present  work. 

We  will  use  a  body-fitted  grid  for  numerical  computation.  This  is  done  by 
applying  the  following  mapping  to  the  physical  domain.  Figure  13  illustrates  the 

mapping  of  a  curved  domain  into  a  rectangular  domain.  Let  y  =  ya[x^  and  y  =  yh(x)  be 

the  bottom  and  top  boundary  of  a  curved  domain  as  shown.  The  mapping  function  to  be 
used  is, 

’  n= y.(!‘)-yXxiy^°~n^+n'y" {x>~'v,|,t*)] 

The  end  points  r\Q  and  rje  are  free  parameters  that  one  may  choose  to  ensure  that  the 
mapped  domain  is  rectangular. 

(xe,yb(xe)) 


ri 


Figure  13.  Mapping  of  a  curved  domain  into  a  rectangular  domain. 


13 


The  mapping  will  transform  the  nozzle  in  the  x-y  plane  into  a  rectangular  domain 
in  the  r\  plane.  So  the  computational  domain  g-p-z  is  a  rectangular  box.  Because 
the  two  side  plates  of  the  nozzle  are  flat  plates,  the  mean  flow  is  two-dimensional.  But 
with  a  3  dimensional  field  of  random  entropy  blobs  as  input,  the  indirect  combustion 
noise  generation  process  is  three  dimensional.  Figure  14  shows  the  mean  flow 
computational  domain  in  the  £-7]  space.  The  height  of  the  nozzle  is  ho.  £  =  0  is  the 
location  of  the  nozzle  throat. 


Figure  14.  The  computational  domain  in  the  plane. 


The  governing  equations  are  the  Euler  equations  and  the  energy  equation.  We  use 
T,  p,  u  and  v  as  dependent  variables.  These  equations  are, 


dT  dT  dT 

—  +  u — + v — + 
dt  dx  dy 


du  dv 
dx  dy 


du  du  du  T  dp 

— +u — + v — + - — 

dt  dx  dy  yp  dx 


0 


dv  dv  dv  T  dp 

— +u — +  v — + - — 

dt  dx  dy  yp  dy 


0 


dp  dp  dp 
-Yt+uTx+vTy+yp 


du  dv 
dx  dy 


For  computing  the  mean  flow,  the  7-point  stencil  Dispersion-Relation-Preserving  (DRP) 
scheme  is  used.  The  DRP  scheme  is  a  highly  accurate  large  stencil  high-order  scheme.  To 
start  the  solution,  the  one  dimensional  solution  is  employed  as  the  starting  condition.  The 
solution  is  time  marched  until  a  steady  state  solution  is  achieved. 
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Figure  15.  The  physical  domain  of  the  upper  half  of  a  military  styled  nozzle. 

Imposing  correct  boundary  conditions  is  extremely  important  for  acoustics 
problems.  For  the  nozzle  under  consideration,  the  inflow  enters  the  computational 
domain  (see  figure  15)  from  the  left  boundary  and  exits  the  domain  through  the  right 
boundary.  Since  the  nozzle  is  symmetric  with  respect  to  the  center  x-z  plane,  it  is 
sufficient  to  compute  the  flow  in  the  upper  half  of  the  nozzle.  With  this  understanding, 
the  following  boundary  conditions  are  enforced. 

a.  Left  inflow  boundary  -  a  hard  boundary  condition  with  T,  p,  u  and  v  -  0  is  imposed. 
Because  all  the  variables  are  prescribed,  outgoing  disturbances  cannot  leave  from  this 
boundary.  They  will  all  be  reflected  back  into  the  computational  domain  and 
eventually  exit  from  the  outflow  boundary  on  the  right  side. 

b.  At  the  top  of  the  computational  domain,  the  wall  boundary  condition  is  imposed. 

c.  At  the  bottom  of  the  computational  domain  shown  in  figure  15,  a  symmetric  boundary 
condition  is  enforced. 

d.  No  boundary  condition  is  needed  at  the  right  or  exit  boundary  (see  figure  15).  This  is 
because  the  flow  at  the  nozzle  exit  is  supersonic.  No  disturbance  can  propagate 
upstream  against  a  supersonic  flow. 

In  our  computational  algorithm,  we  include  shock  capturing  capability.  This  is  because 
the  flow  is  supersonic  downstream  of  the  nozzle  throat.  We  anticipate  the  sharp  nozzle 
throat  would  lead  to  the  formation  of  shocks  inside  the  nozzle. 


Numerical  Results 

Inside  the  nozzle,  the  flow  is  two-dimensional.  In  order  to  display  the  spatial 
distribution  of  the  computed  results,  we  plot  the  pressure  and  the  mean  flow  velocity  u 
along  lines  of  constant  p  .  Figure  16  shows  these  lines  in  the  physical  domain.  Line  0 
is  the  nozzle  symmetric  line.  Line  7  corresponds  to  the  top  wall  of  the  nozzle. 
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Figure  16.  Lines  of  constant  rj  in  the  physical  domain  along  which  computed 

results  are  shown. 


Figure  17.  Pressure  distributions  as  functions  of  x  along  the  various  lines 

shown  in  figure  16. 
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Figure  17  shows  pressure  distributions  in  the  flow  direction  along  the  various  line 
of  figure  16.  Line  7  is  the  pressure  distribution  along  the  upper  wall  of  the  nozzle,  x  =  0 
is  the  location  of  the  nozzle  throat.  This  figure  indicates  that  the  static  pressure  drops  in 
the  flow  direction.  This  pressure  drop  is  the  driving  mechanism  for  a  velocity  increase  in 
the  flow  direction  inside  the  nozzle.  Line  7,  in  particular,  indicates  that  there  is  a  sharp 
drop  in  pressure  until  the  nozzle  throat  is  reached.  Beyond  the  throat,  there  is  a  pressure 
recovery.  If  we  follow  the  pressure  distribution  along  line  4,  the  pressure  recovery  is  very 
abrupt.  This  means  that  a  shock  is  formed.  There  is  a  sharp  pressure  rise  in  all  the 
subsequent  lines,  namely,  lines  3,  2,  1,  and  0.  The  shock  location  propagates  downstream 
and  becomes  closer  and  closer  to  the  nozzle  symmetric  plane.  Figure  18  shows  the 
location  of  the  shock  inside  the  nozzle.  The  shocks  from  the  top  and  bottom  half  of  the 
nozzle  intersect  near  the  nozzle  exit.  They  pass  through  each  other  before  exiting  outside 
the  nozzle. 


Figure  18.  Location  of  shocks  (red)  inside  the  military  styled  nozzle. 


Figure  19  shows  the  axial  distribution  of  mean  flow  velocity  u  along  line  7  of 
figure  16.  There  is  a  sharp  drop  in  u  for  line  0  to  line  4.  They  are  the  corresponding 
velocity  jumps  across  the  shock. 
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Figure  19.  Axial  distribution  of  mean  flow  velocity  u  along  various  lines  of  figure  16. 


Figure  20.  Distribution  of  pressure  contours  inside  the  military  styled  nozzle. 

Figure  20  are  pressure  contours  inside  the  nozzle.  Near  the  nozzle  exit  the 
pressure  contours  merge  into  a  sharp  discontinuity.  This  is  just  another  indication  of  the 
presence  of  a  shock  and  a  transmitted  shock  near  the  exit  of  the  nozzle. 
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It  is  know  that  small  disturbances  such  as  vorticity  waves  will  generate  sound 
when  passing  through  a  shock  wave.  We  expect  the  same  is  true  when  entropy  blobs  are 
convected  through  an  oblique  shock.  This  mechanism  of  sound  generation  has  not  been 
studied  before.  Our  present  investigation  may  help  to  shed  light  on  the  effectiveness  of 
this  shock-entropy  wave  interaction  mechanism  of  noise  generation. 


IV.  Future  research 

Currently  we  are  in  the  process  of  developing  a  parallel  multi-core  computer  code 
for  computing  indirect  combustion  noise  generation  caused  by  the  passage  of  a  stochastic 
entropy  wave  field  through  the  mean  flow  inside  a  military  styled  nozzle.  We  expect  the 
code  development  to  be  completed  soon.  We  will  then  perform  validation  tests.  Once  this 
step  is  completed,  we  will  perform  a  parametric  study  to  determine  the  level,  the 
spectrum  and  other  features  of  the  indirect  combustion  noise  field  that  is  generated  inside 
the  nozzle  of  a  military  aircraft. 
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